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ABSTR  ACT 


A  stable  nuaerlcal  procedure  for  solving  the  coupled* 
nonlinear  equations  of  electron  and  ion  conservation  and 
electrical  potential  (Poisson's  equation) ,  is  described.  A 
tvo-disensional  nodel  with  periodic  active  sites  on  a  flat 
plate  is  utilized  to  obtain  both  qualitative  and  quantita¬ 
tive  results  which  clearly  illustrate  the  self-generating 
sheath  and  anbipolar  regions  adjacent  to  a  non-enitting 
electrode.  The  active  sites  are  characterized  as  voltage 
sources  and  by  electron  densities  depressed  froe  both  the 
non-active  wall  and  the  free-streaa  values.  Various  cases 
of  species  density  at  the  nonactive  wall  are  evaluated. 
Recoabination/ionization  plays  an  important  role  in  estab¬ 
lishing  the  boundary  layer  behavior  in  the  anbipolar  region 
and  in  the  dimensionality  of  the  problen  fornulation. 
Application  is  for  nitrogen  gas  at  one  anagat  in  a  pulsed 
discharge  of  about  50  nicro-seconds  in  duration. 

Results  clearly  deaonstrate  the  boundary  layer  nature  of 
the  species  density  at  the  electrode.  A  sheath  length  of  38 
aicrons  about  a  35-volt  active  site  and  5-20  nicrons  along 
the  non-active  portion  of  the  electrode  is  established. 

Joule  heating  is  deternined  to  be  not  iaportant. 
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I.  lH2£222£Xiafi 


A.  BACKGROUND 

Losses  at  the  electrode  boundary  regions  play  an  impor¬ 
tant  role  in  high  energy  density  aolecular  gas  flow  lasers. 
Of  particular  interest  to  this  study  is  the  Carbon-Dioxide 
glow-discharge  plasaa ,  of  which  Nitrogen  is  the  principal 
energy  bearing  constituent. 

The  nain  energy  loss  aechanisas  in  a  collision  dominated 
plasaa  are  found  in  the  sheath  and  aabipolar  regions.  The 
sheath  voltage  loss  occurs  directly  as  a  result  of  Debye 
shielding,  which  creates  a  space  charge  layer  abutting  an 
electrode  or  wall.  The  aabipolar  losses  are  associated  with 
aabipolar  diffusion,  relating  to  the  charged  species  density 
being  less  than  the  free-streaa  or  eguilibriua  density. 

This  loss  exists  even  is  the  absence  of  current  flow. 

An  initial  preaise  of  probe  theories  that  takes  account 
of  collisions  between  charged  particles,  atoas  and  molec¬ 
ules,  is  that  the  plasaa  is  gua si-neutral  up  to  a  certain 
distance  from  the  probe  and  that  the  electric  fields  are 
localized  within  a  narrow  sheath  region  adjacent  to  the  sur¬ 
face  of  the  probe.  The  problem  of  making  the  transition 
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from  the  quasi-neutral  solution  to  the  sheath  solution  was 
usually  resolved  through  matching  techniques.  The  descrip¬ 
tion  of  a  collision  dominated  sheath  together  with  the  qua¬ 
si-neutral  ambipolar  region  presents  an  imposing  set  of  cou¬ 
pled  nonlinear  partial  differential  equations,  for  which 
there  exists  no  known  exact  solution.  This  study  reports  a 
numerical  method  that  describes  the  sheath  and  ambipolar 
regions  in  a  self-consistent  fashion. 

In  the  hydrodynamic  approximation  (i.e.,  collision  domi¬ 
nated),  the  dimensions  of  the  space  charge  region  "^SH  (the 
sheath  length)  is  larger  than  the  mean  free  path  "/  "  of  the 
species  particles.  Additionally  the  sheath  characteristic 
length  is  very  much  smaller  than  the  characteristic  length 
of  the  plasma  boundary  layer.  This  puts  our  study  in  the 
collisionally  thin  sheath  region,  [3ef.  1]. 

However,  the  plasma  is  not  considered  to  be  in  thermal 
equilibrium.  In  fact  for  plasmas  in  thermal  nonequilibrium, 
the  electron  temperature  is  governed  by  the  electron  energy 
equation,  which  includes  the  effects  of  the  electron  energy- 
conduction  and  electric  field  as  well  as  collisional  energy 
exchange  between  electrons  and  heavy-gas  particles.  The 
consideration  of  negligible  energy  exchange  during  colli¬ 
sions  (thermally  frozen)  does  not  generally  result  in  a 
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constant  electron  temperature  (Tc) ;  particularly  in  the 
sheath  region,  even  if  other  properties  such  as  the  ion 
temperature  (T{)  are  constant.  Therefore  the  electron 
energy  equation  is  necessary  even  in  the  case  of  the  ther- 
mally  frozen  situation. 

An  ionizer/sustainer  discharge  pumped  laser  requires 
that  relatively  high  E/n  values  be  distributed  as  uniformly 
as  possible  throughout  the  bulk  of  the  gas  medium.  The 
electrode  regions  however  represent  a  substantially  higher 
and  more  inhomogeneous  E/n  which  may  impose  an  instability 
on  the  remainder  of  the  discharge.  In  the  region  of  the 
electrode,  the  E/n  is  significantly  higher  than  the  undis¬ 
turbed  plasma.  In  this  region  the  elevated  electron  temper¬ 
ature,  which  is  a  direct  function  of  the  sustainer  field, 
will  increase  ionization,  which  could  lead  to  arcing  and 
breakdown. 

In  the  so-called  vacuum  arc  [Ref.  2],  anode  spots  appear 
under  certain  operating  conditions.  These  anode  spots  are 
of  practical  importance,  having  been  also  observed  in  high 
density  plasmas  [Ref.  3].  They  are  typically  micron  or 
sub-micron  sized  protrusions  on  the  electrode  surface,  which 
act  to  "focus'*  the  electric  field  and  current  stream  lines. 
Primary  consideration  here  is  with  the  anode. 
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B.  REVIEW  Of  PREVIOUS  BORE 
1.  l^isjcs 

Host  work  on  sheath  phenomena  is  eabodied  in  "probe" 
theory  investigations,  where  the  effects  of  the  electrode 
(probe)  locally  disturb  the  quiescent  plasma.  Such  work  is 
relevant  to  this  study  since  the  anode  is  essentially  a 
heavily  biased  probe  in  contact  with  the  plasna,  which  is 
quiescent  within  the  sheath  region. 

Analyses  of  electrostatic  probes  in  collision  doni- 
nated  plasmas  were  originally  done  by  assuaing  that  a  quasi¬ 
neutral  diffusion  controllgfl  region  extended  to  within  one 
aean  free  path  of  the  probe  where  it  was  matched  to  the  edge 
of  the  free  fall-sheath.  Ho  provision  was  aade  for  a  tran¬ 
sition  region  between  the  two  regiaes  (Ref.  4]. 

Shultz  and  Brown  [Ref.  5]  studied  the  transition 
froa  a  collisionless  to  collision  dominated  sheath  through 
charged  species  collection  with  a  spherical  probe.  Postu¬ 
lating  different  physical  models  for  cases  ranging  from  a 
collisionless  sheath  to  a  sheath  with  many  collisions, 
se ai-eapirical  formulas  were  obtained  in  each  instance. 
Several  other  authors  attempted  matching  the  solution  in  the 
collision  dominated  region  directly  to  that  of  the  free-fall 
sheath  using  variational  principles. 


The  previously  mentioned  techniques  can  he  criticised, 
because  in  each  of  these  aethods  something  was  forced  to 
fit.  In  other  words,  at  least  one  derivative  was  disconti¬ 
nuous  at  sone  artificially  imposed  boundary,  thereby  raising 
doubt  as  to  the  validity  of  the  results. 

The  first  systematic  analysis  of  probes  in  collision 
dominated  plasmas  were  carried  out  by  Cohen  [Ref.  6]  and  Su 
and  Lam  [Ref.  7],  who  assumed  the  continuum  equations  valid 
throughout  the  plasma,  including  to  the  probe  surface,  is  a 
consequence  of  these  analyses,  both  the  sheath  and  the  qua¬ 
sineutral  regions  appeared  naturally  from  the  diffusion 
equations  in  the  limit  as  the  ratio  of  probe  size  to  Debye 
length  approached  infinity.  The  continuum  model  removed  the 
concern  about  forcing  sets  of  equations  to  match  at  the  edge 
of  an  arbitrarily  assumed  sheath.  These  authors  restricted 
the  problem  to  an  isothermal  plasma,  but  as  will  be  shown 
here,  the  electron  temperature  may  vary  significantly  in  the 
neighborhood  of  the  probe  for  most  discharges  of  interest. 

Badbill  [Ref.  8]  quasilinearized  the  continuum  equa¬ 
tions,  and  solved  with  numerical  integration,  obtaining 
results  that  agreed  with  those  of  Cohen  where  comparison 


could  be  made 


In  the  previous  analyses,  not  only  was  the  isother- 

■al  plassa  assumed,  but  the  plasaa  was  so  slightly  ionized 

that  only  charge* neutral  interactions  needed  to  be  taken 

into  account.  However,  for  highly  nonequilibriun  plasaas, 

charge-charge  collisions  play  an  iaportant  role  [Refs.  9  and 

10].  The  effects  of  collisions  between  charged  particles 

aanifest  theaselves  as  voluae  ionization,  volume  recoabina- 

tion  and  the  dependence  of  the  transport  coefficients  on 

charged  particle  density  and  temperature  and,  therefore  on 

position.  Barad  [Ref.  11]  cites  a  specific  ezaaple, 

"If  the  electron?  have  a  teaperature  of  2500  ®K,  then  the 
ionization  fraction  equal  to  about  10-*  will  result  in 
charge-charge  collisions  occurring  one- tenth  as  frequently 
as  cnarqe-neutral  collisions.  It  is  shown...  that  for 
such  a  degree  of  charge-charge  interaction,  significant 
effects  are  felt. " 

Chung,  Talbot  and  Touryan  [Ref.  1],  in  an  extensive 
review  of  probe  studies,  state  that  no  general  solution  is 
available  for  deteraining  charge  density  and  species  teaper¬ 
ature  for  probes  saall  in  conpadson  to  the  hydrodynaaic 
boundary  layer  thickness. 

Dolson  [Bef.  12]  investigated  the  nature  and  extent 
of  the  voltage  drops  in  the  vicinity  of  aagnetohydrodynaaic 
(HBD)  non-eaitting  electrodes,  in  particular  the  losses 
attributable  to  the  sheath.  Onder  his  conditions  the  nonex¬ 
istence  of  a  one-diaensional  sheath  solution  is  shown,  and  a 
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computer  aodel  with  two-dimensional,  periodic  active  sites 
representing  a  flat  plate  electrode  is  developed.  With  this 
model  the  effects  of  a  magnetic  field  and  Joule  heating  are 
studied  and  the  results  are  compared  with  experimental 
observations  reported  by  Argryoupolos  (1973) ,  and  Son ju  and 
Teno  (1974),  [Refs.  13,  14). 

An  extension  of  Poison's  work  is  reported  in  1980 
[Ref.  15],  for  application  to  ionizar/sustainer  discharges. 

The  role  of  electron  pressure  is  introduced  along  with  net 

\ 

ionization  and  recombination,  which  proves  necessary  in 
establishing  truely  boundary-layer  behavior  for  the  charged 
species  profiles.  Additionally,  a  modified  Nevton-Raphson 
computational  procedure,  which  enables  much  faster  solution 
convergence  is  formulated.  That  study  demonstrates  solu¬ 
tions  with  non-active  wall  densities  depressed  from  the  free 
stream  density,  a  result  Dolson  did  not  accomplish.  Results 
were  for  10  volts  of  electrode  potential  fall,  an  increase 
from  5.0  volts  reported  by  Dolson. 

2.  saaeataUaaal  asttais 

The  fundamental  difficulty  in  the  development  of 
computational  models  for  the  prediction  of  plasma  flow  phe¬ 
nomena  lies  in  the  coupling  of  equations  governing  both  the 
gasdynamic  flow  and  the  electromagnetic  (EH)  fields.  Since 
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the  flow  fields  can  affect  the  EH  fields  (and  vice-versa), 
the  nonlinear  partial  differential  equations  aust  be  solved 
siaultaneously  and  the  foraulation  of  the  problea  becoaes 
coaplicated.  In  the  present  study  the  flow  properties  are 
decoupled  froa  the  the  charged  species  conservation  equation 
and  Poisson's  Equation. 

The  approaches  used  in  the  solution  of  plasaa  flow 
phenoaena  are  quickly  reviewed  below  so  as  to  point  out  the 
background  of  the  the  approach  in  this  study.  This  problea 
is  typical  to  those  dealing  with  plasaa-boundary  layer  phe¬ 
noaena.  Approaches  to  solving  the  general  plasaa  flow  prob- 
leas  have  consisted  of  the  following; 

(a)  One  diaensional  "influence  coefficients." 

(b)  One-diaensional  closed  fora  analytic  solutions  to 
ordinary  differential  equations. 

(c)  siailarity  transforaation  and  closed  fora  solutions 
to  a  aodified  differential  equation  in  one-  or  two- 
diaensions. 

(d)  One-,  two-,  or  three- diaensional  nuaerical  solution 
to  the  uncoupled  steady  plasaa  flow  phenoaena 
equations. 

(e)  One-,  two-,  or  three-diaensional  nuaerical  solution 
to  the  coupled  steady  plasaa  flow  phenoaena  equations 
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using  either  the  aagnetic  induction  eguation  or  the 
elliptic  eguation  governing  the  electric  fields. 

(f)  One-,  two-,  or  three-dinensional  transient  numerical 
solution  to  the  coupled  plasna  flow  phenoaena 
equations  using  either  the  aagnetic  induction  or  the 
elliptic  eguation. 

Methods  (a)  through  (c)  are  exact  but  are  restricted 
to  a  narrow  range  of  probleas.  Approaches  (d)  and  (e)  pro¬ 
vide  rapid  convergence  to  a  steady  state  solution  but  do  not 
consider  the  transient  behavior  between  the  gasdynaaic  and 
electroaagnetic  fields,  or  often  the  strong  coupling  between 
the  two  sets  of  fields.  Method  (f)  is  aost  desirable  since 
it  treats  a  wide  range  of  probleas,  but  suffers  froa  large 
coaputational  tines  [  Ref .  16].  The  approach  of  this  study 
is  characterized  as  a  aethod  (e)  approach. 

C.  OBJECTIVES  OP  THE  PRESENT  STUDY 

The  objective  of  this  study  is  to  deaOnstrate  a  nuaeri- 
cal  procedure  which  can  nodel  the  physical  behavior  of  a 
conducting  fluid,  as  would  be  generated  in  the  after-glow  of 
an  electron  bean  discharge  laser,  under  the  influence  of 
electric  forces.  Calculations  are  perforaed  in  two  diaen- 
sions,  with  a  flat  plate  non-eaitting  electrode  (anode) 
which  has  periodic  active  sites.  These  active  sites 
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facilitate  current  constriction  required  to  satisfy  charge 
continuity,  Ohe's  lav,  and  Poisson* s  equation.  The  physical 
sodel  is  idealized  in  Fig.  1. 

The  controlling  equations  are  solved  for  the  potential, 
the  charged  species  density  distributions,  electron  tempera¬ 
ture  and  current  in  the  vicinity  of  the  electrode.  The 
resulting  computer  programs  generate  sheath  and  aabipolar 
regions  in  a  self-consistent  fashion,  using  the  saae  set  of 
equations  throughout  the  field.  The  size  of  the  sheath  and 
the  voltage  drop  attributable  to  its  shielding  effects  are 
investigated. 

This  work  does  aot  consider  those  phenoaena  that  lead  to 
the  arc  discharge,  which  results  from  thermal  instabilities. 
The  constriction  at  the  anode  is  purely  that  required  to 
satisfy  continuity. 
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ii.  paoBiea  fobapiahoh 


A.  GOVEBNING  EQUATIONS 

The  governing  equations  that  describe  the  electrode 
region  portion  of  the  plasaa,  consist  of  Poisson's  equation, 
species  conservation  and  electron  energy  equation,  in  addi- 
tion  to  the  overall  continuity,  aoaentua  and  energy  equa- 
tions  of  the  plasaa.  These  overall  plasaa  equations  are 
essentially  those  of  the  background  gas  since  the  plasaa  is 
weakly  ionized  (less  than  0.01  %) . 

Ne  have  [Bef.  1]; 

Poisson  equation... 

electron  conservation . . . 

tt- Vne  -  v  (S * 

ion  conservation... 

(4l7p;  ♦  * 
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In  gas  discharges,  the  electron  energy  distribution  is 
known  to  be  non-Naxwellian.  In  order  to  allow  for  non-Hax- 
wellian  distributions  and  still  retain  use  of  a  teaperature 
concept,  an  empirical  fit  to  the  calculated  relation  of 
electron  temperature  as  a  function  of  E/n  is  used  in  the 
representation  of  the  electron  teaperature  (see  Sec.  III. A) . 
The  concept  of  a  "two-temperature"  plasaa  is  preserved  but 
no  restrictions  are  made  on  the  energy  distribution  function 
of  the  electrons.  In  addition  to  the  assumption  of  a  func¬ 
tional  dependence  of  the  electron  teaperature  upon  E/n  only, 
there  is  no  representation  of  an  external  magnetic  field 
and,  in  particular  the  internal  or  self-generated  aagnetic 
fields  are  not  considered.  The  overall  gas  teaperature  is 
assumed  constant  and  in  eguilibriua  with  the  plasaa  ions. 

The  bulk  energy  equation  is  assumed  to  have  little  bearing 
on  this  study,  since  there  is  negligible  gas  heating  as  dis¬ 
cussed  in  section  III. A.  This  effectively  uncouples  overall 
energy  consideration  froa  the  problem. 

The  electrodes  serve  the  purpose  of  sustaining  the  S/n 
to  puap  the  laser  aediua.  No  account  is  aade  of  an  external 
E-beaa  or  other  source  of  priaary  ionization.  Siaply,  the 
aediua  aay  be  considered  to  exist  in  an  afterglow  state. 
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The  production  tera  (i.e.,  the  right  hand  side)  in  Eqns. 
(2)  and  (3)  has  been  described  utilizing  aeasured  and  theor¬ 
etical  results  of  the  rate  coefficients  for  Nitrogen  ioniza¬ 
tion  and  two-,  and  three-body  recoabination.  The  production 
tera  is  expressed  as; 

ne  =  ^ne-*nenf  ayient.  n0  (4) 

where  and  "d*"  are  respectively  the  ionization  and 

recoabination  rate  coefficients. 

B.  BOONDABI  CONDITIONS 

In  the  custoaary  boundary- value  problea  in  fluid  aechan- 
ics,  a  two-diaensional  flat  plate  is  assuaed.  A  siailar 
two-dimensional  Cartesian  description  is  used  in  the  present 
work,  which  has  been  found  to  be  the  ainiaua  suitable 
description  for  the  problea  posed,  [Refs.  3,12,15].  No 
one-diaensional  solution  exists  for  the  frozen  flow  of 
charges  through  a  collisional  sheath,  since  a  Cartesian  geo- 
aetry  in  one  diaension  is  devoid  of  the  necessary  geoaetric 
decrease  of  currant  density  away  froa  the  electrode,  as 
there  would  be  in  a  cylindrical  or  spherical  geoaetry. 

Convection  is  considered  negligible,  as  can  be  deaons- 
trated  by  coaparing  field  velocity  with  the  drift  velocity. 
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Electros  and  ion  densities  are  determined  at  the  plasma 
edge,  knowing  the  state  of  the  aediaa.  A  current  balance 
can  be  used  to  establish  the  species  densities  resident  at 
the  active  portion  of  the  anode,  An  electrode  voltage  drop 
is  deternined  froa  the  properties  of  the  aediua;  e.g.,  the 
effective  ionization  energies  of  the  gas  in  particular. 

Since  current  constricts,  there  nust  be  regions  of  the  flat 
plate  that  receive  no  current.  The  inactive  portion  of  the 
flat  plate  aerits  attention,  especially  in  the  regions  adja¬ 
cent  to  the  active  anode  site.  Species  concentration  along 
the  inactive  portion  of  the  wall  can  vary,  depending  on  the 
reactivity  (catalytic  properties)  of  the  wall.  There  is  a 
balance  that  can  be  established  between  the  voltage  varia¬ 
tion  noraal  to  the  wall,  and  the  species  concentration  at 
the  wall  together  with  the  species  variation  noraal  to  the 
wall.  A  zero  concentration  of  electrons  and/or  ions  at  the 
wall  is  an  acceptable  boundary-condition;  which  is  to  say 
that  the  wall  acts  as  a  perfect  sink  for  the  charged  spe¬ 
cies.  A  detailed  kinetic-theory  analysis  shows  that  there 
is  in  fact  a  lower  bound  to  the  electron  and  ion  concentra¬ 
tions  at  the  wall  [ Bef .  17].  In  any  case,  there  is  no  cur¬ 
rent  into  or  out  of  the  inactive  wall,  (see  Fig.  1) .  This 


is  discussed  later  in  Sect.  IV 


The  remaining  portion  of  the  boundary  conditions  refer 
to  the  free  stream  medium,  and  are  easily  established. 

C.  PARAMETERS  OP  THE  DISCHARGE  HEDIOH 

The  application  of  this  study  is  for  molecular  Nitrogen 
at  one  amagat,  where  the  gas  temperature  is  273  °K.  Other 
pertinent  parameters  follow; 

Electric  Field:  1.5  to  15  KV/cm. 

Electron  Density:  10** 7  to  10**»  m”3. 

Electrode  Voltage  Drop:  35  volts. 

The  electron  energy  relation  is,  as  mentioned,  a  curve  fit 
from  experimental  results  [Ref.  18],  carried  out  by  R.R. 
Comptoa  and  D.J.  Sutton  (1952).  Since  cross-section  data 
for  electrons  in  Nitrogen  discharges  are  available  (Ref. 
19],  the  electron  temperature  and  diffusion  coefficient  are 
readily  determined,  and  representable  as  a  function  of  E/n. 

The  ion  temperature  is  assumed  to  follow  that  of  the 
neutral  gas.  An  electrode  potential  fall  of  about  35  volts 
may  be  assumed  since  it  is  the  effective  ionization  energy 
for  Nitrogen  gas  and,  therefore,  the  anode  fall  [Ref.  20]. 

The  cathode  is  more  complicated  to  describe  because  of 
the  requirement  of  surface  electron  emission.  The  cathode 
fall  in  cold  cathodes  is  of  the  order  of  250  volts,  and  is 
considerably  less  in  thermionic  cathodes.  A  proper 
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III.  PROBLEM  §otas£OM 


A.  PROBLEM  SIMPLIFICATION  (PHYSICS) 

There  are  a  number  of  s imp lit ic at ions  that  are  incorpo¬ 
rated  in  the  problem  description.  Namely; 

•  Steady  state  conditions  (within  the  pulse  time). 

•  No  magnetic  fields. 

•  Negligible  convection. 

•  Negligible  Joule  heating.  =  constant. 

However,  specific  account  is  made  of; 

•  Ionization/recombination. 

•  Electron  pressure. 

These  are  discussed  in  part  in  this  section. 

1 .  go  &  vest  jo  a 

The  ambipolar  region  is  the  transition  region  fcjm 
the  sheath  to  the  undisturbed  plasma  and  can  perhaps  span 
the  boundary  layer.  This  fact  requires  some  consideration 
of  the  presence  of  convective  effects. 

From  the  ion  species  conservation  equation  (Eq.  3). 
for  example,  we  know  that  in  the  ambipolar  region  the  last 
term  of  the  left  hand  side  is  small  compared  to  electric 
conduction.  Now  we  need  only  to  compare  the  convection  with 

d 
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the  conduction,  or  simply  the  fluid  velocity  with  the  drift 
velocity. 

.ic 

vs^,B  «  0.1  to  1.0  a/sec 

where  the  Reynolds  nuaber  (Re)  is  104  to  10*  and  the  free 
stream  velocity  is  10  to  100  a/s.  Now  the  drift  velocity 
aay  be  approximated  by 

Vj  as  10  to  100  a/s 

for  E  *  10s  to  10*  V/a  and  the  mobility  /l  3  10-4 
a*/s/v olt.  Clearly,  in  the  presence  of  a  relatively  strong 
interelectrode  field,  the  contribution  of  ion  convection  due 
to  a  cross  flow  should  be  negligible.  The  geometrical 
orientation  of  the  discharge  with  respect  to  the  priaary 
flow  can  thus  becone  an  important  consideration. 

2.  Daasiiz  aauodatx  Uki  Eataatiaa  Sins 

The  charged-species  density  profiles  change  in  a 
more  complicated  manner  than  either  the  voltage  potential  or 
the  electric  field.  The  magnitude  of  the  ion  and  electron 
densities  aay  change  appreciably  within  the  sheath;  aore- 
over,  fractional  analysis  is  risky  because  inflections  are 
present  in  the  profiles,  A  discussion  pertaining  to  the 
stability  of  the  aabipolar  region  and  its  boundary  layer 
nature  is  presented  in  detail  by  Riblarz,  et  al. ,  [Ref.  15]. 
However,  Appendix  1  shows  that  the  characteristic  length  of 
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the  aabipolar  region  is  of  the  order  of  the  characteristic 
sheath  length.  Letting  L=5x10“*  asters  and  the  drift  veloc¬ 
ity  as  before,  v^  *  10  to  100  a/s,  a  characteristic  tiae  for 
the  foraation  of  the  density  boadary  layer  would  be  0.5  to 
5.0  aicroseconds.  It  is  reasonable  that  the  aabipolar  dif¬ 
fusion  to  the  waLl  can  establish  itself  well  within  the  50 
aicrosecond  pulse  period  of  the  electrodes. 

3.  siaaif Isaacs  QZ  Jaals  isatiaa 

For  the  present  problea,  joule  heating  can  be  shown 
not  to  be  significant.  In  the  vicinity  of  the  active  site, 
the  average  energy  deposition  rate  per  unit  voluae  is  shown 
to  be  appro xiaately  4x10*  watts/a3.  The  teaperature  of  the 
Nitrogen  gas  at  one  aaagat  would  be  raised  approxiuately  1°K 
in  2.6x10-*  sec.  Allowing  a  50x10-*  sec  pulse  period  would 
raise  the  teaperature  of  that  region  by  approxiuately  20°K, 
or  7*. 

Assuae  that  the  heat  generated  at  the  anode  spot  is 
carried  away  by  theraal  diffussion  (conduction).  Then  for 
Nitrogen  gas  (k*0.0267  N/a-°K,  /  *1.138  kg/u3,  and  €^*1043 
<J/kg-°K)  the  theraal  diffusivity  is; 

fl(s  •pjjp  *  2.25x10“*  a*/sec 

A  typical  diffussion  tiae  would  be  (for,  *4.4x10-*  a) ; 

1 0.86  aillisec 
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The  "flush  tiee"  is  the  tiee  required  for  clearing 
the  lasing  cavity  of  residae  gases  in  preparation  for  the 
next  pulse,  and  can  be  esti sated  as  the  gas  flow  rate 
divided  by  the  flow  length.  Assuming  v  =  100  a/sec  and  L  = 
10  ca,  the  flush  tiee  would  be; 

FT  =*  0.1/100  *  1  aillisec 

other  factors  considered,  the  inter-pulse  period  (tiee  bet¬ 
ween  pulses)  would  be  larger,  but  even  in  a  one  Billisecond 
period,  there  is  just  enough  tiae  to  cool  the  anode  spot. 

The  analysis  so  far  has  assumed  a  flow  velocity  at 
the  surface  to  be  the  sane  as  the  bulk  flow  in  the  laser 
cavity.  However  a  boundary  layer  exists —  probably  turbu¬ 
lent,  with  a  laeinar  sub-layer  adjacent  to  the  electrode 
for; 

8C  3  3  (0.1)  (100)/(  15.53x10-*)  *  6.44x10* 

The  velocity  boundary  layer  thickness  would  be  2.55  ailliae- 
ters.  Assuming  a  1/7-th  power  profile,  the  velocity  at  a 

sheath  length  froa  the  electrode  would  be: 

Vi 

v  *  1^(A,/f)  3  100(4.4x10-8/2.55x10-*)  V7  =  55  e/sec 
A  characteristic  local  flush  tiae  would  be  approximately  2 
ailliseconds,  and  auch  higher  for  a  laainar  velocity  pro¬ 
file.  These  flush  tines  can  be  quite  large  as  the  require¬ 
ment  to  cool  the  Joule  heated  gas  near  the  anode  spot  is 


increased.  The  design  flush  tiae  based  on  other  factors  aay 
infact  be  insufficient  to  clear  enough  of  the  Joule  heated 
gas  froa  the  vicinity  of  the  anode  spot,  enabling  a  build  up 
of  heated  gas,  which  could  lead  to  local  enhanced  ionization 
of  the  gas. 

4.  Electron  Temperature 

The  electron  teaperature  (Te )  is  assuaed  to  be  a 
linear  function  of  E/n  as  presented  in  Pig.  2,  for  the 
region  of  interest,  naaely  E/a  *  0.5  to  5.0x10-*®  ?-a*.  Por 
the  range  of  E/a  considered,  it  is  adeguate  to  take  the 
electron  teaperature  as; 

T€  =  [A*log(E/n)  ♦  b]*1^  (5) 

where; 

k  *  12.  1 


B  3  38.5 

B/n  is  normalized  to  10** 0  V-a* 

The  diffusion  coefficients  were  deterained  froa 
Brown  (Hef.  18].  The  drift  velocities  were  obtained  froa 
plotted  results,  and  since  /E,  and  D^v^kT^/e,  (Ein¬ 

stein's  relation) ,  the  diffusion  coefficients  are; 

0^  *  0.0665  a* /sec 
Dfc.  *  8. 0x10-*  a*/sec 
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B.  BOB KING  EQUATIONS 


Since  convection  is  negligible,  and  assuming  that  the 
sustainer  operates  in  the  afterglow  of  an  electron  beaa,  the 
governing  equations,  nanely  Egs.  1,2  and  3,  can  be  written 


as  follows 


*  c,  (ne-  n£) 


(6) 


-7*[  Vn^ ♦  /&  -  n eW/*]  =crne  (7) 

C  AA  a aA  a 

-7«[7n*  ♦  £n;V0]  »  c,nc  (8) 

where; 

n3/n^  $-  e^/hT^ 

6s  re/Tt«*  <*  =  T«^/T0 

c,  =  e^/lcT^  ct=  Aj2nt^/De  c3a  C2.D*/Dc 

The  equations  nay  now  be  applied  to  flat  plate  fornula- 
tion  as  depicted  in  Fig.  3.  It  is  convenient  to  drop  the 
"•a"  for  simplicity,  but  it  is  to  be  understood  that  all 
variables  have  been  suitably  non-diaensionalized. 


Be  have  Poisson's  equation; 


electron  conservation  equation; 


J  0  |5x  d* 


5*j  5*j 


1 


3x* 


Sfl 


(10) 


*€. 


ion  conservation; 


c,fA;(w«-n  0 


C,fle  (11) 


C.  MOMERICAL  METHODOLOGY 

The  nonlinear  terns  creat  problems  in  the  coaputer  anal 
ysis  of  the  electron  and  ion  species  conservation  equations 
(Bqs.  10,11).  The  Jacobi  aethod  includes  all  nonlinear 
terms  on  the  "right  hand  side";  i.e.,  external  to  the  coef¬ 
ficient  aatrix.  Convergence  to  a  solution  is  possible  if 
these  non-linear  terns  change  slowly  over  each  iteration. 

It  was  found  by  Dolson  [Ref.  12]  that  the  Jacobi  aethod  was 
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in  fact  unstable  for  the  present  set  of  equations  and 
conditions.  As  a  consequence  of  the  failure  of  the  Jacobi 
aethod,  a  quasi-Jacobi  eethod  vas  iapleaented,  in  which  an 
estinate  for  each  of  the  solution  variables  {$  ,  na,  n^)  was 
coeputed.  When  the  product  of  two  variables  is  encountered* 
one  variable  is  treated  as  a  constant  coefficient  for  each 
iteration.  This  aeans  that  the  non-linear  teras  are 
retained  in  the  coefficient  aatrix.  The  "constant"  coeffi¬ 
cients  are  updated  after  every  iteration,  thus  altering  the 
coefficient  aatrix.  The  conventional  Jacobi  net hod  was 
found  to  converge  only  for  low  values  of  voltage  potential 
at  the  electrode,  whereas  the  quasi-Jacobi  procedure  pro¬ 
vided  converged  solutions  for  approxiaately  5  volts.  A  dis¬ 
cussion  of  these  aethods  is  provided  in  Appendix  C  as  back¬ 
ground  for  the  developaent  of  the  aodified  Newton-Raphson 
(MIR)  aethod,  which  is  used  in  the  solution  to  Eqs.  1 
through  4. 

The  nonlinear  coupled  partial  differential  equations 
(POE),  Eqs.  9,10  and  11  are  linearized  to  provide  a  set  of 
linear  algebraic  equations,  which  aust  hold  throughout  the 
doaain  described  as  follows.  The  doaain  as  illustrated  in 
Pig.  3,  is  divided  up  into  equal  sided  sub-areas  whose  ver¬ 
tices  are  the  coaputational  nodes  at  which  the  linearized 
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algebraic  equations  aust  hold.  The  wall  is  represented  by 
seven  (7)  nodes,  with  the  leftmost  node  being  the  location 
of  the  active  site.  The  nyN  direction  extends  from  this 
active  site  along  the  wall.  The  "xw  dimension  perpendicular 
to  the  wall  is  represented  by  25  nodes  starting  from  each  of 
the  wall  nodes.  The  solution  to  the  PDB's  is  represented  by 
the  vector  z(x,y)  ,  or  on  the  computational  mesh  as  z(i,j). 

The  boundary  conditions  have  been  discussed  previously 
and  are  presented  in  Fig.  3.  The  boundary  condition  that 
there  is  no  net  current  into  the  nonactive  wall  (^=0)  is 
satisfied  by  the  electron  contribution  to  current  only. 

This  is  justified  since  it  can  be  shown  that  the  electron 
contribution  is  approximately  three  orders  of  magnitude 
greater  than  the  ion  contribution  to  the  total  current. 
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17.  SSSSiXS  MD  COHSL05JON5 


A.  GENERAL  DISCUSSION 

Computer  solutions  ace  reported  here  for  three  cases 
presented  in  Table  I.  The  active  site  voltage  vas  35  volts 
in  each  of  the  three  cases.  Included  in  this  section  are 


TABLE  I:  Coaputer  Solution  Conditions 


CASE 

ne  (vail) 

nt*  (vail) 

n*  (anode)  n^  (anode) 

I 

1.0 

O 

• 

o 

1 

1 

J-natched  0.0 

II 

"float" 

o 

• 

o 

not  J-natched  0.0 

III 

"float" 

"float" 

not  J-aatched  0.0 

also  soae  pertinent  results  obtained  at  10  volts  to  illus¬ 
trate  the  effects  the  net  production  tera  has  cn  the  conver¬ 
gence  to  a  boundary  layer  behavior  for  the  species  density 
and  free  streaa  electric  field. 

The  solutions  are  depicted  as  obligue  presentations  of 
the  potential,  species  density,  space  charge,  net  produc¬ 
tion,  electric  field  and  Joule  heating  distributions  over 
the  "r-y"  coordinates  defined  in  the  problea  foraolation. 
Also  included  are  tvo-diaensional  plots  of  various  terns 
along  a  “cut"  or  line  perpendicular  to  the  electrode  vail. 
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starting  at  the  active  site,  with  soae  starting  at  the  vail 
along  the  inter-nodal  line  of  symmetry  (see  Pig.  3) .  The 
solutions  are  organized  as  Case  I, II  and  III  results,  mill 
be  presented  as  Pigs.  12,13  and  14  respectively  at  the  end 
of  this  section.  The  discussion  of  these  results  proceeds 
with  additional  figures  which  are  interspersed  in  the  text 
appropriately. 

Table  I  presents  various  electron  and  ion  boundary  con¬ 
ditions  along  the  wall  soae  of  which  are  "floated.”  By  this 
tera  is  meant  that  upon  the  completion  of  each  nuaerical 
iteration,  the  electron  (and/or  ion)  density  at  the  nonac¬ 
tive  portion  of  the  wall  is  re-established  by  assigning  it 
the  value  of  a  near  neighbor  away  from  the  wall.  This  pro¬ 
cess  is  illustrated  in  Pig.  4. 
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Figure  4 

Description  of  "floating  the  Nall  Boundary  Condition. 
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The  usual  elliptic  boundary  value  problem  has  well 


possed  boundary  conditions  from  a  priori  information.  How¬ 
ever,  the  nature  of  the  problem  in  this  study  includes  very 
little  information  as  to  the  correct  species  wall  density. 

It  seems  reasonable  to  set  these  densities  to  zero  along  the 
inactive  portion  of  the  electrode  as  previously  mentioned  in 
Sec.  II. B.  In  this  study  the  wall  boundary  conditions  are 
"lossened"  to  facilitate  "naturally  occuring"  and  unpresum¬ 
ing  species  density  profiles.  The  nature  of  the  problem  in 
effect,  produces  a  self-consistent  boundary  condition. 

B.  EFFECTS  OF  NET  PRODUCTION 

Earlier  work  [Ref.  IS]  reported  limited  results  for  a  10 
volt  active  node  site.  Solutions  for  cases  with  and  without 
the  inclusion  of  the  net  production  term  have  profound 
influences  on  the  appearance  of  the  charged  species  density 
profiles.  Figures  5,6  and  7  illustrate  the  before  and  after 
effects  of  including  the  net  production  term.  In  both  cases 
the  wall  boundary  values  for  the  species  density  was  ini¬ 
tialized  at  the  free-streaa  value  (i.e.,  1.0),  then  allowed 
to  "float"  as  described  previously.  The  addition  of  the  net 
production  term  brings  about  the  boundary  layer  behavior  of 
the  densities  of  the  electrons  and  ions  is  illustrated  in 
Fig.  7.  Efforts  to  attain  converged  solutions  with  the  net 
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j  Figure  7 

Influence  of  the  net  production  tern  on  the 
density  profiles  along  a  line  perpendicular  to  the  vail, 
I  extending  froe  the  active  node  site,  (10  volts) . 


C.  DISCUSSION  OF  RESULTS 


JfeUaaa.  Ealsaiiai 

The  voltage  potential  of  the  active  node  site  ("the 
anode1*)  in  each  of  the  cases  presented  was*  35  volts.  Fig¬ 
ures  12-a,  13-a  and  14-a  will  clearly  illustrate  sharp  vol¬ 
tage  spikes  at  the  anode,  and  the  "Laplacian"  nature  of  the 
field  far  fron  the  anode.  The  potential  along  the  inactive 
vail  varied  froa  12  volts  at  the  farthest  position  froa  the 
anode  to  17  volts  just  next  to  the  anode.  Sinilar  profiles 
of  the  electric  potential  along  the  vail  for  the  other  cases 
are  presented  in  Fig.  8.  Figures  12-i,  13-i  and  14-i  vill 
shov  the  voltage  potential  variation  along  a  line  noraal  to 
the  vail  extending  froa  the  anode  site. 

Figure  9  illustrates  a  tvo-diaensional  cut  of  the 
voltage  potential  along  the  inter-nodal  syaaetry  line.  Here 
it  can  be  seen  that  the  variation  froa  the  vail  is  shallow. 
It  is  iaportant  to  point  out  that  the  electron  teaperature 
is  directly  related  to  the  gradient  of  the  voltage  potential 
through  Eg.  5.  This  fact  has  certain  conseguences  which  are 
discussed  in  the  next  section. 
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Pigure  8 

Voltage  potential  along  the  wall  boundary  (electrode) 
for  CASES  I,  II,  and  III. 
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Figure  9 

Voltage  potential  variation  normal  to  the  Mali 
along  tne  inter-nodal  symmetry  line,  (CASE  I) 


2.  Electron  Energy 


The  electron  energy  will  be  presented  in  series  -b 
of  Figs.  12,13  and  14,  as  oblique  perspectives.  The  anode 
site  tenperature  for  Cases  I,  II  and  III  are  2.1,  2.1  and 
1.9  e?  respectively.  The  description  of  the  electron  temp¬ 
erature  is  that  of  a  linear  function  of  the  logarithm  of  the 
electric  field.  Any  portion  of  the  potential  field  that  is 
"level**  (grad  0  >0)  will  necessarily  generate  a  low  or  even 
negative  (see  Fig.  14-b)  electron  teeperature.  According  to 
Fig.  2  and  the  discussion  in  Sec  II. A. 4,  a  negative  tempera¬ 
ture  would  occur  below  an  E/n  of  0.04x10“*°  V-cm*.  This 
problem  can  be  removed  by  setting  a  minimum  value  to  the 
non-dimensional  temperature  field,  nominally  0.13,  corres¬ 
ponding  to  the  minimum  at  the  low  end  of  E/n  in  Fig.  2. 

The  parameters  of  the  discharge  (Sec.  I.C)  require 
an  electric  field  of  2.7x10*  V/m  (or,  log£B]s5.43)  at  the 
free  stream  boundary,  which  corresponds  to  an  electron  temp¬ 
erature  of  0.905  eV  (<?*1.).  The  temperature  fields  in  all 
cases  presented  do  not  quite  meet  the  free  stream  boundary 
requirements.  The  electric  fields  are  insufficient  to  ren¬ 
der  proper  results  to  the  electron  temperature  profiles. 

Thus  the  free  stream  temperature  boundary  requirements  are 
"forced"  by  slightly  inflating  the  influence  of  the  electric 
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field.  The  effect  of  forcing  the  teaperature  field  in  this 
way  can  also  be  looked  at  as  "sliding"  the  linear 
representation  of  the  teaperature  curve  to  higher  values  at 
constant  E/n,  (see  Fig.  2).  Figure  10  illustrates  a  case 
where  a  low  electric  field#  here  0.8  kT/ca#  results  in  a  low 
electron  teaperature  projecting  toward  the  free  streaa 
(0*0.62).  The  low  teaperature  is  adjusted  to  aeet  the 
desired  free  streaa  value  (0*1.0)  by  aultiplying  the  elec¬ 
tric  field  by  a  factor  of  3.40.  Pigure  10  also  illustrates 
an  alternate  aethod  of  aeeting  the  free  streaa  conditions  by 
truncating  the  variation  of  teaperature  with  a  ainiaua 
value#  naaely  the  free  streaa  boundary  condition. 

Both  of  the  "forcing"  techniques  are  used#  but  for 
different  reasons.  &s  long  as  the  adjustaent  factor  does 
not  exceed  an  order  of  aagnitude  ange  in  the  electric 
field#  then  it  is  not  objectionable  in  light  of  the  siapli- 
fied  representation  of  the  electron  teaperature.  The 
adjustaent  reguireaents  of  the  three  cases  reported  in  this 
study  are  between  3.0  and  5.0#  with  Case  I  using  a  factor  of 
10  to  see  the  effects  on  the  Joule  heating.  In  Case  I  the 
non-diaensional  teaperature  is  artificially  high  at  1.4  eV 


at  the  free  streaa 


NON- CM  ME  NSl  ON  AL  TEMPERATURE  (6) 


DISTANCE  FROM  THE  ACTIVE  MODE  SITE 

*  5 

(multiplies  of  X  *  4.4  10  m) 


Figure  10 


Electron  temperature  variation  noraal  to  the  vail,  along  the 
nodal  syaaetry  line;  illustrating  tvo  methods  of  adjustment 
to  "force"  the  free  stream  temperature  conditions. 
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The  electric  field  near  the  inter-nodal  symmetry  points 
of  the  wall  ate  characteristically  2  orders  of  aagnitude 
less  than  that  of  the  free  streaa  conditions.  The  resulting 
very  low,  even  negative  teaperatures  in  Case  III,  Fig.  14-b, 
is  "patched"  using  the  truncation  aethod  previosly  dis¬ 
cussed.  whether  these  low  teaperatures  and  the  associated 
steep  derivatives  are  physically  reasonable,  merits  further 
investigation. 

The  temperature  fields  generated  and  presented  as 
Figs.  12,  13,  14  series  -b  are  the  final  result  depictions. 
However,  the  temperature  fields  used  in  the  solution  process 
were  adjusted  and  truncated  as  previously  discussed  to  main¬ 
tain  solution  convergence.  An  example  of  a  "corrected" 
temperature  field  is  presented  in  Fig.  14-c. 

3.  ascitis.  ziald  sZaaia  t&aiiaa 

The  electric  field  and  Joule  heating  distributions 
are  generated  from  the  primary  voltage  and  electron  density 
solutions,  much  in  the  Banner  as  the  temperature  distribu¬ 
tion  was  determined.  In  fact,  the  Electric  field  and  the 
joule  heating  distributions  are  siailar  in  appearance  to  the 
temperature  distributions.  These  distributions  play  no 
direct  part  in  the  solution  scheme  and  thus  being  ancillary 
in  nature  were  not  "forced"  to  meet  any  a  priori  conditions. 


The  Joale  heating  term  can  be  written  as; 

J.E  *  [  De  n#ae)c3^  /^*  ].[^  v*  -tfne  ]*V^ 
and  thus  is  a  rather  straight-forward  function  of  4,  i^, 
and  a; .  The  ion  contribution  is  shown  to  be  negligible, 
since  Dfc»  Dc- .  The  Joule  heat  at  the  anode  site  is  of  order 
10*  and  is  100  tiaes  that  of  the  free  streaa  values.  Joule 
heating  effects  are  negligible  in  the  overall  problem,  as 
discussed  in  Sec.  II. A. 3,  obviating  the  reguirement  for  an 
energy  eguation. 

4.  SR&siss.  Psasjti  ttisttifeaiiaa 

The  solutions  for  the  charged  species  density  dis¬ 
tributions  are  presented  as  oblique  perspectives  in  Figs. 
12,13  and  14,  series  -e  (electron  density);  series  -f  (ion 
density)  ;  and  series  -g  (space  charge  density)  .  Addition¬ 
ally  the  variation  of  these  distributions  noraal  to  the  wall 
along  the  active  node  symmetry  line  are  presented  in  the 
series  -j  and  -k  figures.  The  oblique  perspectives  of  the 
densities  are  reversed  for  clarity. 

The  electron  and  ion  densities  at  the  nonactive  wall 
are  a  subject  of  interest  in  the  literature,  to  which  little 
information  has  been  gathered  on  &  firm  basis.  The  results 
of  this  study  suggest  a  slightly  positive  space  charge  at 
the  inactive  wall.  Scrutiny  of  Case  III  results  reveals,  a 


2.2x10**  «• 3  positive  space  charge  number  density,  (or  3.5 
milli-coulombs) .  Recall  that  Case  III  involved  "floating" 
both  of  the  charged  species  densities  at  the  vail.  In  fact, 
close  inspection  of  Case  I  and  II  reveals  a  positive  space 
of  similar  magnitude  but  displaced  slightly  from  the  wall. 
The  wall  in  these  cases  was  set  with  a  negative  space  charge 
as  the  boundary  condition.  The  positive  space  charge  along 
the  wall  is  most  predominant  in  the  vicinity  of  the  inter- 
nodal  symmetry  line,  and  whether  the  low  temperature  values 
in  this  region  have  significance,  remains  to  be  determined. 

The  active  node  site  (anode)  electron  density  in 
Case  I  (2. 4x10* 7  n-3)  was  determined  so  as  to  match  the  free 
stream  current  to  the  current  constricting  at  the  anode.  In 
the  results  of  Cases  II  and  III  the  anode  electron  densities 
are  chosen  so  as  to  provide  a  smooth  variation  extending 
outward  from  the  active  site. 

5.  Ionization/Recoab ina tio  n 

The  net  ionization  tern  is  utilized  in  all  cases; 
and  was  found  necessary  for  (1)  establishment  of  the  charged 
particle  density  distributions,  and  (2)  solution  convergence 
for  sheath  lengths  spanning  more  than  one  computational  node 
spacing.  The  net  production  tern  is  depicted  in  Figs.  12,13 
and  14  series  -h.  There  is  net  ionization  throughout  the 


solution  domain.  However  as  Case  III  demonstrates,  there  is 
a  factor  of  5  orders  of  aagnitude  greater  ionization  within 
the  sheath  region  surrounding  the  anode  site  than  at  the 
free-stream  edge.  Cases  I  and  II  show  ionization  along  the 
nonactive  wall,  siaply  because  the  ion  density  is  zero,  eli- 
ainating  significant  recoabination,  (see  Eg.  4). 

6.  cacEgfli  ass2ti£ii2B 

Figure  11  is  typical  of  the  constricting  nature  of 
the  current  streaa  lines.  The  potential,  electron  density 
and  the  electron  temperature  are  used  to  deteraine  the  "x-y" 
components  of  the  current  density.  Using  a  “nearest  neigh¬ 
bor”  influence  principal,  a  local  derivative  can  be  deter¬ 
mined;  froa  which  a  streaaline  can  be  "aarched-out”  in  saall 
increments.  The  profiles  of  all  three  cases  are  siailar, 
the  fact  being  that  the  current  aatching  discussed  previ¬ 
ously  will  not  affect  the  current  streaaline  profiles.  Fig¬ 
ure  11  clearly  demonstrates  the  satisfaction  of  the  current 
constriction  reguireaent  necessary  to  perait  solutions  to 
the  systea  of  eguations. 

The  current  into  the  anode  site  in  Case  I  is  matched 
to  the  current  entering  at  the  free-streaa  boundary.  This 
specifies  the  electron  density  at  the  anode  site.  However, 
in  Cases  II  and  III  the  current  at  the  anode  was  only  three 


52 


tines  the  free-stream  current.  The  current  aatching  problem 


does  not  deserve  serious  complaint,  since  the  problem 
description  is  simplified  to  tvo-liaensions  and  the  current 
constriction  is  only  required  to  satisfy  continuity. 


Figure  11 

Two-dimensional  depiction  of  the  current  stream  line 
constriction  toward  the  active  node  site. 
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Figure  12-e 

Oblique  perspective  of 
Electron  Density  (CASE  I) 
(reversed  perspective) 


Figure  12-g 

Oblique  perspective  of 
Space  Charge  (CASE  I) 


Figure  12-f 

Obligue  perspective  of 
Ion  Density  (CASE  I) 
(reversed  perspective) 


Pigure  12-h 

Oblique  perspective  of 
Net  production  (CASE  I 
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Figure  12-i 

Presentation  of  the  Voltage  Potential 
(CASE  I)  along  a  line  nornal  to  the  »»1 
extending  froe  the  active  node  site. 
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Figure 

Presentation  of 
(CASE  I)  along  a  line 
extending  troa  the 
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Figure  12-1 

Presentation  of  the  Electric  Field 
(CASE  I)  along  a  line  noreal  to  the  wall, 
extending  froe  the  active  node  site, 
(logaritha  scaled) 
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Figure  12-s 

Presentation  of  the  Joule  Heating 
(CASE  I)  along  a  line  norsal  to  the  wall, 
extending  ffroa  the  active  node  site, 
(logarithm  scaled) 
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Figure  1  3-a 


Oblique  perspective  of 
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Figure  1 3-c 


Ooligue  perspective  of 
log  Electric  Pield  (CASE  II) 


Figure  13-b‘ 

Oblique  perspective  of 
Teaperature  field  (CASE  II) 


Figure  13-d 

Oblique  perspective  of 
log  Joule  Heat  (CASE  II) 


Figure  13-e  Figure  13-f 

obligua  perspective  of  Oblique  perspective  of 

Electron  Density  (CASE  II)  Ion  Density  (CASE  II) 

(reversed  perspective)  (reversed  perspective) 


Figure  13-g 

Oblique  perspective  of 
space  Charge  (CASE  II) 


Figure  13-h 

Oblique  perspective  of 
Net  Production  (CASE  II) 
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Figure  13-j 

presentation  of  Species  Density 
(CASE  IIV  along  a  line  nornal  to  the  wall, 
extending  froa  the  active  node  site. 
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Figure  13-k 

Presentation  of  the  negative  Space  Charge  Density 
(CASE  III  along  a  line  nomal  to  the  wall, 
extending  froa  the  active  node  site. 
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Pigure  13-1 

Presentation  of  the  Electric  Fi*l* 
(CASE  III  along  a  line  noraal  to  the  wall* 
*  extending  froa  the  active  node  site, 
(logaritha  scaled) 
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figure  13-e 

Presentation  of  the  Joule  Heating 
(CASE  II)  along  a  line  nornal  to  the  v 
extending  from  the  active  node  site 
(logaritha  scaled) 
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Figure  13-n 

presentation  of  the  Species  Densities 
(CASE  III  on  a  line  normal  to  the  mall, 
extending  along  the  Inter-nodal  symmetry  line. 
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Figure  14-a  Figure  14-b 


Figure  14-c 


Figure  14-d 


Obligue  perspective  of  adjusted 
Temperature  Field  (CASE  III) 


obligue  perspective  of 
Joule  Heat  TCASE  TTT1 


Figure  14-e 

EiSsys'Disfuniiviu* 

(reversed  perspective! 


Figure  14-f 

Obligue  perspective  of 
Ion  density  (CXSE  III) 
(reversed  perspective) 


Figure  14-g 


Oblique  perspective  of 
Space  Charge*^  (CASE  III) 


Figure  14-h 


Oblique  pe 
Met  Product 


s 


71 


0.00 


0.  10 


0.  15 


0.05 

DISTANCE  (MM) 


Presea tatio 
(CASS  III1  al 
extending 


Pigure  14-i 
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Figure  14-j 

Presentation  of  Species  Density 
(CASE  III)  along  a  line  noraal  to  the  wall, 
extending  froa  the  active  node  site. 
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Pig  are  14-lc 


Presentation  of  the  negative  Space  Charge  Density 
(CASE  III)  along  a  line  noraal  to  the  vail, 
extending  from  the  active  node  site. 
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Pigure  14-e 

Presentation  of  the  Joule  Heating 
(CASE  IIIV  along  a  line  noreal  to  the  wall, 
extending  froe  the  active  node  site, 
(logarithe  scaled) 
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Figure  14-n 

Presentation  of  the  Species  Densities 
(CASE  III)  on  a  line  nornal  to  the  wall, 
extending  along  the  inter-nodal  synaetry  line. 


V.  CONCLUSIONS  &ND  RECOMMENDATIONS 


A.  REVIEW  OF  ASSUMPTIONS 

In  the  developnent  of  this  work  several  assuaptions  were 
aade  which  the  solution  results  later  verified.  The  assump¬ 
tion  of  negligible  convection  effects  is  straight  forward. 
The  charged  particle  density  boundary  layer  evolution  tiae 
turns  out  to  be  snail  in  relation  to  the  pulse  period;  thus 
enabling  consideration  of  the  problea  as  a  steady-state  aft¬ 
er-glow  problen.  The  fact  that  the  Joule  heating  of  the  gas 
is  shown  negligible,  frees  the  foraulation  froa  an  elaborate 
energy  description.  Since  the  electron  tenperature  is  a 
siaple  function  of  E/n,  the  problen  can  be  further  sinpli- 
fied.  The  solution  proceeds  with  noainal  values  of  the  par- 
aaeters  of  the  discharge  as  discussed  in  Sec.  II. B. 

Finally,  the  iaportance  placed  on  the  net  production  tera  is 
verified  by  the  profound  changes  in  the  results  including 
the  fact  that  "n^"  is  necessary  for  the  establishaent  of  the 
charged-species  boundary- layer  behaviour.  The  analysis  of 
App.  3  shows  that  an  assuaption  of  "frozen  flow"  would  in 


fact  be  preaature 


B.  PHYSICAL  CONC LU SIGNS 

The  present  sod el  of  the  sheath  and  aabipolar  region 
evolves  fro a  the  assuoptions  of  steady  state,  unifora  ion 
teaperature  and  constant  diffusion  coefficients.  A  suaaary 
of  iaportant  results  is  presented  in  Table  II. 


TABLE  II:  SUMMARY  OF  IMPORTANT  RESULTS 


SUMMARIZED 

RESULTS 

CASE 

I 

CASE 

III 

(units) 

sheath  length 

3.8 

3.5 

3.7 

10-*  m 

aabipolar  length 

4.5 

3.5 

3.7 

n 

r 

1 

o 

E^(free-  stream) 

7.9 

3.9 

6.3 

10  ♦  V/m 

3,  (anode) 

1.9 

2.5 

2.2 

10*  V/a 

The  results  show  that  a  sheath  develops  as  anticipated 
froa  theory  and  the  proposed  assumptions.  The  sheath  char¬ 
acteristic  length  is  consistent  with  that  predicted  by  Debye 
shielding.  Debye  shielding  of  a  flat  plate  predicts  that 
the  sheath  length  is  4.4x10"*  aeters.  For  a  point  (or 

spherical)  shielding  the  sheath  length  is  3.1x10-*  aeters. 
The  present  study  involves  the  shielding  of  a  point  voltage 
source  situated  on  a  surface  charged  plate.  Since  both  two- 
and  three-diaensional  effects  are  inherent  in  the  formula¬ 
tion,  the  resulting  sheath  lengths  are  quite  reasonable. 
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As  Table  II  shows,  the  size  of  the  aabipolar  region  is  of 
the  saae  order  as  the  sheath  size.  Here  the  ambipolar 
length  is  defined  as  the  length  extending  froa  the  edge  of 
the  sheath  to  the  point  at  which  99X  of  the  free-streaa  den¬ 
sity  is  recovered.  Pinally,  the  sheath  can  be  self-generat¬ 
ing  froa  a  consistent  set  of  egaations  without  an  overall 
energy  equation. 

The  resulting  electric  fields  at  the  free-streaa  boun¬ 
dary  do  not  guite  attain  the  desired  2.7x10s  7/a  presented 
in  the  problea  formulation.  These  results  fall  short  by  the 
saae  factors  required  for  the  adjustaent  of  the  electron 
teaperature  field  to  Beet  the  stated  free-streaa  boundary 
reguirenents.  It  aust  be  pointed  out  that  the  electric 
field  is  the  required  boundary  condition.  However,  this 
boundary  condition  could  be  aade  to  yield  no  reasonable 
results,  and  so  the  chosen  boundary  condition  is  6  -  0. 

The  resulting  electric  fields  at  the  active  site  in  each 
of  Cases  I,  II,  or  III  was  of  two  orders  of  aagnitude  higher 
than  that  of  the  free-streaa.  These  electric  fields  are 
below  that  needed  for  breakdown  of  Nitrogen  at  one  atmo¬ 
sphere,  which  is  deterained  as,  [Ref.  18]; 

2*  3  (2/p)  p  »  (350  V/ca-torr)  x(760  torr)  *  2.7x10'  7/a 
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The  requirement  of  strictly  Batched  currents  produces 
steep  derivatives  and  inflections  of  the  density  profiles  in 
the  vicinity  of  the  anode.  I  noticeable  "spike"  is  evident 
in  Pig.  12- j  for  example.  Case  II  and  III  anode  electron 
densities  were  established  on  more  aesthetic  grounds,  where 
saoother  appearing  density  profiles  were  generated.  In 
these  cases  the  anode  current  exceeded  the  free  streaa  cur¬ 
rent  by  only  a  factor  of  3. 

The  nonactive  regions  of  the  electrode  (wall)  far  fron 
the  active  site  reveal  a  snail  positive  space  charge  den¬ 
sity  aaounting  to  approxiaa tely  231  of  the  free-streaa  den¬ 
sity.  The  negative  space  charged  sheath  is  46%  to  56%  of 
the  free-streaa  density. 

Previous  work  [Ref.  12]  of  a  similar  problem  reported  no 
reasonable  results  for  the  case  of  a  catalytic  wall 
(iijSD;  sO).  The  results  presented  in  this  study  deaonstrate 
wall  conditions  of  -  0.  and  <  1 .  Here,  a  noticeable 
positive  space  charge  coexists  with  the  large  negative  space 
charged  sheath,  without  large  inflections  in  either  the  vol¬ 
tage  potential  or  densities. 

The  Case  III  study  has  both  electron  and  ion  wall  condi¬ 
tions  floated.  The  shape  of  these  density  distributions  is 
quite  symmetrical  about  the  anode  site  in  all  directions. 
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One  needs  little  imagination  to  visualize  that  the  density 
profiles  can  be  "wrapped"  around  the  anode  completely,  in 
effect  treating  the  anode  as  a  whisker  or  protrusion  into 
the  plasma.  However,  the  imagination  must  be  tempered  by 
the  fact  that  it  can  only  be  a  two-dimensional  protrusion 
(e.g. ,  a  "blade")  • 

The  effect  of  the  net  production  term  has  been  previ¬ 
ously  discussed  as  one  of  the  most  important  results  of  this 
study.  The  term  is  necessary  for  inclusion  in  the  problem 
in  order  (1)  to  attain  wall  electron  density  values  of  less 
than  the  free-stream,  (2)  to  produce  the  boundary-layer 
behavior  of  the  species  densities, and  (3)  to  obtain  con¬ 
verged  solutions  with  voltages  greater  than  10  volts  or  com¬ 
putational  node  spacings  of  any  order  less  than  that  of  the 
sheath  length. 

C.  NUMERICAL  CONCLUSIONS 

The  Modified  Newton-Rap hson  (MNR)  method,  presented  to 
solve  the  particular  set  of  coupled  non-linear  partial  dif¬ 
ferential  equations  in  this  study,  is  a  general  method  that 
can  be  applied  to  more  sophisticated  systems.  The  defi¬ 
ciency  of  the  quasi- Jacobian  method  (discussed  in  Sec.  III.C 
and  App.  C.)  in  dealing  with  non-linear  terms  is  circum¬ 
vented  by  an  alternative  linearizing  scheme.  Convergence  is 


82 


rapid;  however,  the  starting  solution  oust  not  exceed  10 
volts,  oust  have  zero  net  production,  and  oust  have  a  cons¬ 
tant  electron  tenperature  during  the  first  2-3  iterations. 
After  the  solutions  take  their  initial  form  at  tan  volts, 
the  potential  field  is  "amplified''  and  reinstated  as  a  new 
starting  solution.  The  final  stages  of  the  converging  solu¬ 
tion  are  carried  out  under  the  full  influences  of  (1)  float¬ 
ing  wall  densities,  (2)  variable  electron  temperature,  (3) 
net  production  and  (4)  an  active  node  site  electron  density 
being  solved  for  through  current  matching. 

The  profiles  of  the  space  charge  and  net  production  dis¬ 
tributions  demonstrate  erratic  behavior  at  the  first  few 
computational  nodes  from  the  active  site.  This  is  believed 
to  be  caused  by  the  linearization  of  the  governing  non-li¬ 
near  equations  and  has  the  appearance  of  truncation  error 
when  observed  on  an  iteration  to  iteration  basis.  However, 
the  remaining  regions  of  the  profiles  are  smooth. 

0.  RECOMMENDATIONS 

This  work  succeeded  in  solving  the  problem  in  a  two-di¬ 
mensional  environment.  This  fact  necessarily  limits  the 
quantitative  results  since  the  current  density  decreases  in 
only  two  dimensions.  Kith  the  computational  power  of  recent 
computers,  it  should  prove  worthwhile  to  make  the  next 
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advance  as  a  three-dimensional  problem  rather  than  elaborate 
in  the  two- dimensional  case,  A  three-dimensional  scheme 


would  render  more  realistic  results  of  a  quantitative  nature 
because  active  spots  are  more  likely  to  exist  than  active 
lines  and  because  3-0  geometries  require  less  voltage  for  a 
given  current. 

Specific  recommendations  other  than  the  above  include; 


•  Account  for  the  ion  contribution  to  the  current 
boundary  condition  into  the  wall. 

•  Investigate  the  physical  implications  of  the  low 
electron  temperature  near  tne  wall  far  from  the 
anode  site. 

•  Investigate  the  physical  implications  of  the  ele£- 
tron  density  "spikes"  prevalent  near  the  anode  site. 

•  Hake  the  wall  temperature  a  valid  boundary  condi¬ 
tion. 

•  Ose  a  free-stream  boundary  condition  of 

*2.7x10*  V/m,  rather  than  zero  potential.  This 
would  introduce  desired  results  for  pumping  lasers. 

•  Include  a  more  sophisticated  electron  temperature 
model,  or  verify  the  present  one. 

•  Increase  the  computational  mesh  density. 

•  Study  the  general  application  limits  of  the  HNB 
method. 

•  Investigate  "n  "  in  the  role  of  making  one-dimen¬ 
sional  solutions  possible  (i.e. ,  alter  the  reacting 
rates) . 

•  Investigate  whether  the  low  electron  temperatures 
and  associated  steep  derivatives  located  along  the 
wall  between  active  sites  have  any  physical  implica¬ 
tions,  or  result  merely  from  inadequate  problem  for¬ 
mulation. 

•  Consider  the  electron-beam  ionization  in  the  solu¬ 
tion  to  produce  more  realistic  effects. 
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APPENDIX  A 


CHAB ACTERI STIC  PARAMETERS 

A.  SHEATH  LENGTH 

The  extent  of  the  sheath  is  one  of  the  sost  important 
charcteristic  lengths  in  this  study.  It  is  within  the 
sheath  that  sost  of  the  potential  drop  occurs  for  a  short 
discharge.  Fortunately,  the  sheath  length  can  be  estimated 
rather  easily  froa  Poisson's  equation,  (Eg.  AM 


7V  *  £j(we-ru) 


Let ; 


*  *  "MST  *  '•  "5  ’  "s/nt* 


(AM 


(A2) 


where  *^11'  is  the  characteristic  sheath  length;  the  sub¬ 
script  "o"  indicates  the  value  at  the  electrode  and  "  "  the 

value  at  the  undisturbed  plasaa.  After  gathering  the  diaea- 
sional  quantities  on  one  side,  Bq.  A2  becoaes; 


The  characteristic  sheath  length  H  Xf"  •  defined  in  this 
fashion,  is  used  to  non-diaensionalize  the  "del"  operator  in 
all  directions.  That  is,  independent  of  the  diaensionality 
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of  the  problea,  the  sheath  is  effective  over  the  character¬ 
istic  length  jast  defined.  Now  an  estiaate  of  the  size  of 
the  sheath  can  aade  by  setting  the  grouping  of  diaensional 
quantities  in  Eg.  A3  to  be  of  the  order  one,  naaely; 


As 

L*««* 


■f- 


where  "  ?I|J  is  the  Debye  length  in  the  undisturbed  plasaa. 

Anticipated  values  are  now  established  for  the  sheath 
length.  Take  for  ezaaple  a  35  volt  anode  potential,  and  a 
charge  species  density  of  10 43  to  104*  a-3.  These  paraae- 
ters  define  a  sheath  length  of  1.4x10~*  to  1.4x10"*  aeters. 

As  can  be  seen,  a  ten-fold  increase  in  voltage  would 
change  "  by  a  factor  of  3,  rendering  the  nondiaensionali- 
zation  of  the  problea  relatively  insensitive  to  the  chosen 
value  of  the  voltage  potential.  That  is  to  say  the  value  of 
the  sheath  characteristic  length  as  defined  aay  be  consid¬ 
ered  a  resonable  estiaate  for  the  charge  species  densities 
exaaplif iad. 

The  iaportant  iaplication  here  is  that  extreae  care  aust 
be  devoted  to  electrode  surface  preparation  to  preclude 
whiskers  or  flaws  of  a  size  coaparable  to  N^N.  These  pro¬ 
trusions  would  produce  a  focusing  of  the  electric  potential 
to  create  extreaely  high  electric  fields,  which  becoae  sites 
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for  initiation  of  breakdown.  Protrusion  growth  has  been 
observed  to  fora  in  a  tiae  on  the  order  of  Billiseconds* 
tlan;  conditions  are  favorable  to  protrusion  growth ,  depen¬ 
dent  on  the  pre-breakdown  and  breakdown  current  profile, 
electrode  coaposition  (especially  in  the  case  of  a  ther- 
aionic  cathode)  ,  surface  purity  and  surface  preparation 
[Bef.  21].  Soae  possibilities  are; 

•  Electrode  plastic  flow  or  aeltinq  in  an  electric 
field  due  to  aicro-particle  boabardaent. 

•  sticking  of  aicroparticles. 

•  Electron  beaa  induced  protrusion  growth. 

B.  DIFPUSION  LENGTH 

If  diffusion  is  a  significant  factor  in  this  study  then 
its  contribution  to  the  current  would  be; 

*  *  D‘*&  - 

or,  since  >  0^,  and  using  Einstein's  relation,  or 
D^4C*kTe/e; 

J  *  =  A  kTe1^ 

Baking  the  equation  unit  free  (through  fractional  analysis) 

-r-r  *  ~  (»5) 

e  Be***  «9x 


97 


If  the  right  side  of  Eg.  AS  is  of  the  order  one,  (the  nondi- 
aensional  density  gradient  is  a  aaxiaua  of  1)  then  "L"  nay 


be  estiaated  as; 

L  a  eDen^/J  sc  3.7x10-*  a 

for  the  values; 

n«  *  10* •  a"3 

*  0.0665  a*/(V*sec) 

=^e»nt«V#  =  J4f+e» nc«nA*  (E/n)  ss  2.88x103  A- a-3 
Since  L,  then  any  role  attributed  to  diffusion 

should  be  visible  within  the  sheath.  The  E/n  within  the 

thin  sheath  is  of  course  greater  than  the  E/n  at  the  free 

streaa  (10*®  7*a*),  which  decreases  the  local  value  for  "L" 

still  further. 

Another  approach  is  to  assuae  a  one-diaensiona 1  solution 
to  the  charged  species  conservation  equations,  which  can 
only  be  done  outside  of  the  sheath,  i.e.,  in  the  aabipolar 
region.  The  ion  conservation  equation  aay  be  written; 

n  *  7«J  ■  eV«[7n^  “  D^n]  (A 6) 

Since  it  can  be  shown  (next  section)  that  the  production 
tera  can  be  approxiaated  as; 

n  atP4*n  -0Cn* 

then  for  the  stationary  condition  of  no  net  production; 

•  ef  N 

Hondiaensionalizing  Eg.  A6,  it  can  be  shown  that  the 
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characteristic  diffusion  length  "LM  is; 

I  *  { D/5Fn  *  2.x10“*  meters 

which  is  of  the  sase  order  of  aagnitude  as  the  characteris¬ 
tic  sheath  length. 

C.  IONIZATION/BECOH BINATIO V  LENGTH 

Analysis  of  the  production  tera  (Eg.  4)  reveals  that  for 
Nitrogen  the  three-body  (electron)  collisional  recoabination 
rate  (#*)  is  saall  in  relation  to  to  the  other  terms,  (see 
App  B.)  .  IT  any  case  the  net  production  tera  is; 

“e  •*}«*-<$»«*  -aJW11#  (i7) 

In  the  steady  state  the  net  production  tern  is  zero,  and 
as  is  done  in  Appendix  B,  the  relation  bee ones; 

s' 

ne  *  (*,♦  n*)*(nr-  n«)  «n€ 

substituting  the  value  of  0. 17x10“* ‘  aVsec,  and  (t?  -  ly)  at 
10**,  we  can  estiaate  the  characteristic  production  tiae  to 
be  Tp  *  5.88x10-*  sec. 

To  establish  a  aeaningful  characteristic  length,  the 
average  theraal  velocity  will  be  used  since  it  represents 
the  principal  transport  velocity  of  the  aediua. 

•  5.88x10-**3.2x10*  *  1 . 88x10“4  a 
The  characteristic  length  for  net  production  is  a  factor 
of  four  larger  than  the  sheath  length  (A$*  4.4x10“*  a).  The 


conclusion  which  can  be  aade  fro*  this  (fractional)  analysis 
is  that  the  net  production  is  significant  within  the  overall 
region  of  interest,  including  the  sheath  itself. 

D.  FBOZEH  FLOS 

Shen  a  process  relaxation  tiee  "7^"  is  long  coapared 
with  the  tiae  that  a  fluid  particle  spends  in  a  given 
region,  the  flow  is  said  to  be  "frozen"  in  this  region.  By 
this  tera  one  neans  that  a  reaction  process  nay  be 
neglected,  as  far  as  influencing  the  flow  is  concerned,  and 
the  species  involved  treated  as  cheaically  inert,  [Bef.  22]. 
In  the  opposite  liait,  where  "7J."  is  very  saall  relative  to 
"4^"  the  characteristic  flow  tiae,  the  flow  is  said  to  be  in 
"eguilibriua. "  For  the  recoabination  rates  given  in  this 
section,  "7J"  *  6.x10"T  and  "7^"  »  4.0x10"*  seconds.  Sere 
flow  tiae  is  defined  as  the  sheath  characteristic  length 
divided  by  the  drift  velocty.  Thus  the  flow  is  not  "froxen" 
as  far  as  cheaical  reaction,  but  we  aust  consider  that  the 
electron  and  ion  species  are  not  in  theraal  egulibriua. 


APPENDIX  3 


DEVELOPMENT  OF  THE  PRODUCTION  TERN 

The  net  production  tern  is  defined  as  the  difference 
between  ionization  and  reconbination  rates.  Whenever  a  neu¬ 
tral  atoa  or  aolecule  loses  one  or  aore  of  its  electrons, 
ionization  has  ocurred.  Recoab ination  is  the  process  wher¬ 
ein  positive  ions  and  electrons  coabine  to  fora  neutral 
atoas  or  aolecules.  It  is  siaply  reverse  ionization.  If 
the  net  production  is  zero  then  ionization  and  reconbination 
are  equal.  The  net  production  tern  is  the  right  side  of 
Eqs.  2  and  3,  the  expression  for  which  is  developed  here. 

This  work  considers  only  three  types  of  yeggi bination 

(1)  Dissociative  -  the  electron  recoabines  with  a  aolecu- 

lar  ion  (AB)  +  ,  and  the  reconbination  energy  goes 
into  dissociating  the  aolecule  and  increasing  the 
kinetic  energy  of  the  resulting  products.  The  pro¬ 
cess  is  described  as; 

e  ♦  (IB)  ♦  — >  A  ♦  B 

(2)  Three- body  (elect sop)  -  the  energy  released  in  the 
reconbination  process  is  carried  off  by  a  third  body 
which  is  an  electron.  The  process  is  described  as; 

e  ♦  k+  ♦  e  A  ♦  e 
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(3)  Thy^e-body  Uia&Hl  “  the  recombination  energy  is 


carried  off  by  a  third  body  which  is  a  heavy  parti¬ 
cle  such  as  an  atos  or  nolecule.  The  process  is 
described  as; 

e  ♦  A+  ♦  B  — ♦  A  ♦  B 

This  study  does  not  consider  other  recoabination  mechan¬ 
ises  such  as  radiative  or  dielectric  recoabination,  since 
they  are  generally  of  low  significance  in  relation  to  the 
just  described  processes.  The  ions  are  liaited  to  singly 
positive  ions.  Ho  account  is  aade  for  interaediate  or 
excited  states  as  this  reguires  a  detailed  balancing  which 
is  beyond  the  scope  of  this  work.  Additionally,  no  account 
is  aade  of  attachaent. 


Bewriting  Bq.  a  for  siaplicity  of  discussion; 

ne  »  if i  ne  -  AJn^a^  •jgiutain0  (B1) 

where  each  "44*  refers  to  the  recoabination  coefficients  for 
dissociation,  three- body  (electron)  and  three-body  (heavy)  , 
respectively. 

The  recoabination  coefficients  depend  upon  the  ionic 
species  involved.  Soae  general  but  representative  values 
[Ref.  22]  follow; 


Hinnov  and  Hirschberg  (  Bef ,  23]  developed  a  relation  for 
three-body  (electron)  recombination  coefficients  as  a  func¬ 
tion  of  electron  temperature  (1^)  and  density,  valid  for 
energies  less  than  0.25  eV. 

1. 09rl0-*oa^T<-*/*  mVsec 

where  the  electron  temperature  is  ®K.  Gurevich  and  Pitaen- 
ski  (Ref.  24]  extended  this  development  to  obtain; 

2. 13x1 0-2® neT€-®/»lnA-  m3/sec 

where; 

lnwi*  1 . 24  x107  (Tq3/ne)>* 

At  1.0  eV  these  coefficients  are  5.6x10~*i  and  1.05X10-1® 
mv/sec  respectively.  A  simplification  of  the  later  expres¬ 
sion  for  is  valid  from  energies  of  0.1  to  10.0  eT  and 
ne=10-»»  a-*; 

Ot-j,  *0. 0  32xT~4*  3 

Brown  [Ref.  18:  p.197]  gives  a  value  for  the  dissocia¬ 
tive  recombination  coefficient; 

C(t  »  1.5x10“12  m*/s 

The  ionization  coefficient  can  be  determined  from  con¬ 
sidering  that  there  is  no  net  production  in  the  free  stream, 
i. e. ,  far  from  the  electrode.  Here; 

0  -*,N*  -4^3  -*jH2n4  (B2) 
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where  "NM  is  the  stationary  state.  Solving  for  and  sub¬ 
stituting  into  Eg  31  gives; 

nft=  n^(N-n'  )  (*i  +  tijn»)  ♦  *1neT$»/2* (H*-i^n4- )  (83) 

where; 

*,  =  1.5x10“ 12  (fro«  Ref.  18) 

tf^n^  1  .Q9x  10“20  (from  Ref.  23) 

=  2 .0x1 0”13  (froa  Ref.  22) 

The  net  production  tern  is  the  relation,  suitably  non- 
diaensionalized,  that  appears  in  the  coaputer  subprogram 
"RTSIDE. " 


APPENDIX  C 


NUMERICAL  SOLUTION  METHOD 


A.  GENERAL  DISCUSSION 


There  are  two  basic  approaches  used  in  solving  a  nonli¬ 
near  systea  of  equations,  such  as  possed  in  this  study.  A 
given  method  can  usually  be  classified  as  one  of  the  ascent 
(or  descent)  or  more  commonly  refered  to  as  gradient  meth¬ 
ods,  or  as  one  of  the  Newton  methods.  The  method  applied  in 
this  study  and  by  Dolson  [Ref.  12]  is  one  of  the  Newton 
methods.  A  brief  review  of  both  methods  will  be  presented 
to  facilitate  development  of  the  particular  method  employed 
in  this  study. 

Let  the  nonlinear  system  of  equations  be  written  as; 
f/  (x|»xJ#****xn)  =  0 


(Cl) 


(xi  »xt*  •  •  •  3  0 

or  in  vector  form  for  convenience; 

f (x)  =0 


(C2) 


In  the  descent  method  the  problem  is  possed  in  such  a 
way  so  as  to  fora  a  function  which  is  zero  at  any  solution 
and  positive  otherwise.  A  typical  such  function  would  be; 
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n 

0  (xi,X2  ,  •  •  .xn)  -^Inl2  <c3> 


so  that  the  problea  is  reduced  to  finding  a  ainiaua  of  0, 
which  also  happens  to  be  the  value  zero.  In  solving  prob¬ 
lems  by  this  method,  an  initial  guess  to  the  solution  must 
be  provided,  call  it  x<l>.  Then  a  direction  d  is  determined 
so  as  to  point  to  a  better  solution,  x<2>.  This  direction 
is  often  chosen  to  be  the  gradient  of  0(x) .  The  operation 
-grad  0(x)  is  the  direction  in  which  0(x)  decreases  most, 
and  is  locally  the  best  direction  in  which  to  reduce  0(x) . 

d  =  -grad  0(x) 


so  that ; 


x<2>  s  x« » >  -  a*grad  0(x<l>) 

where  “a"  is  chosen  to  be  an  appropriate  step  size.  A  prob¬ 
lem  is  encountered  at  local  minima  of  0(x) ,  which  needs  to 
be  avoided.  One  must  also  use  the  step  "a"  in  such  a  manner 
as  to  control  the  acceleration  along  slowly  or  rapidly  vary¬ 
ing  0(x).  It  may  be  that  by  the  nature  of  0(x)  ,  the  gra¬ 
dient  may  have  to  be  approximated  by  finite  differences,  an 
approach  which  offers  a  flexible  recourse  even  when  not 
required. 
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The  Newton  Method  involves  approximating  the  nonlinear 
problem  by  a  suitable  linear  problea  involving  derivatives 
of  the  nonlinear  equations  evaluated  at  x<*>.  The  linear¬ 
ized  problea  is  siaply  a  set  of  linear  algebraic  equations, 
the  solution  to  which  furnishes  an  increaent  that  when  added 
to  the  solution  x<‘>  will  yield  x<*>.  This  Method  is  a  gen¬ 
eralization  to  n-diaensions  of  Newton's  method  for  evaluat¬ 
ing  the  zeros  of  a  function. 

Bach  function  fi  (xj  is  approximated  by  the  initial  two 
terms  of  a  Taylor  series  expansion  about  the  solution  vector 
x<  i>; 

fi  (x)  *  fi(x<*>)  ♦  J  (x<  i>)  •(  x  -  x<‘>]  (C4) 

where  "J  '  is  the  familiar  Jacobian  matrix. 


n 

J-l 


afj(x(,)  ) 


The  right-hand-side  of  equation  C4  is  the  linear  vector 
function  (tangent  hyperplane)  of  x  which  best  approximates 
the  nonlinear  function  f  evaluated  at  x<*>.  Ne  set  this 


right-hand-side  equal  to  zero  and  solve  for  x,  which  is  an 
improvement  from  x<*>.  Me  call  this  improved  solution  vec¬ 
tor  x<*>,  and  from  equation  C4  we  get; 

x<*>  *  tCl)  -  J(X<  »>) -»*f  (X)  (C6) 
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The  Jacobian  need  not  be  known  very  accurately,  and  often  is 
left  unchanged  for  several  iterations  since  its  evaluation 
often  entails  penalizing  computer  tine.  Thus  the  aethod 
(slightly  aodified  as  described  briefly  in  section  III.B) 
used  by  Dolson  [Ref.  12],  is  referred  to  in  aany  texts  as 
the  Hewton  aethod  or  as  the  Hewton-Raphson  aethod  [Refs.  25, 
26]. 

The  Jacobian  aatrix  of  coefficients  for  elliptic  prob- 
leas  is  typically  a  sparse  banded  aatrix,  to  which  there  are 
aany  cob put er  subroutines  available  for  solution. 

B.  BO 01 FI BD  NEiTOH-RAPHSOH  METHOD 

The  particular  solution  technigue  employed  in  this  study 
aerits  a  detailed  discussion,  since  it  is  not  coaaon  in  the 
literature.  Bailey  and  Touryan  [Ref.  27]  use  a  siailar 
aethod  and  provide  soae  description  of  its  developaent.  The 
aethod  is  siailar  to  other  Net  '.on  aethods  however,  the  Jaco- 
bian  aatrix  is  replaced  by  a  aatrix  of  coefficients  devel- 
oped  froa  a  sequence  of  linear  differential  equations  which 
approxiaate  the  systea  of  nonlinear  equations.  This  aethod 
will  be  referred  to  as  a  aodified  Mewton-Raphson  aethod  or 
HMR  aethod  for  short,  and  will  be  presented  in  general  using 
a  specific  exaaple. 


Take  Po_ -son's  equation  given  earlier  and  repeated  here 
for  simplicity; 

V*0»  =,  (n^  -n£)  (C7) 

The  known  solution  variables  are; 

nt<‘>,  and  (>> 

In  general  the  figs^  step  is  linearization  of  the  equa¬ 
tion.  In  this  case  Eq.  C7  is,  linear,  but  it  will  serve  as 
an  exanple.  ie  assuae  that  the  suceeding  solution  variables 
at  iteration-  (2) ,  namely  z<*>#  can  be  approximated  by  the 
solution  z<»>  plus  some  small  change  Az<u. 

Az<*>  a  Z<»>  ♦  Z<»>  (C8) 

where; 

z<*>  =  ty<K>  n .<♦>)*** 

It  remains  to  determine  Az<‘>.  Substitution  of  this 
relationship  into  the  Poisson  Equation  is  the  second  step: 

vV’  r  vV't  rW'~  .  c,  (W-A**') 

rearranged; 

V*A$0,-C, (AWe  -  )  *  -  [tfV*'-  C,  (  «e°-  («> 
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The  right  hand  side  of  Eg.  C9  is  a  function  of  the  known 
solution  vector  *<'».  We  are  applying  this  relationship  at 
each  computational  node,  h  second  order  partial  differen¬ 
tial  eguation  can  be  approx iaated  by  a  second  order  differ¬ 
ence  equation  or  "3-point  formula,"  so  Eg.  C9  will  becoee  a 
system  of  linear  algebraic  relationships  involving  near 
neighboring  values  of  the  solution  vector  z<*>.  The  left 
side  of  Eg.  C9  is  composed  of  a  system  of  coefficients 
expressible  in  terms  of  times  the  unknown  (perturba¬ 

tion)  vector  Az<*>.  So  we  now  have; 

*  -P(Z(D)  (CIO) 

Here  the  matrix  of  coefficients  is  a  sparse  matrix 
with  a  relatively  narrow  band  of  non- zero  elements  about  the 
main  diagonal,  provided  that  the  "natural"  order  of  the 
resulting  linearized  equations  is  preserved.  There  are 
standard  computer  routines  available  which  take  advantage  of 
this  banded  structure,  such  as  LEQT1B,  an  IHSL  routine. 

It  follows  that  since  =  z<*>  - 

z<*>  «  z<i>  -  r  A(z<  »>)  ]-»*P(Z<»>)  (Cl  1) 
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For  the  case  of  a  nonlinear  tera  such  as; 


n« 

there  is  no  additional  difficulty.  When  linearized  this 
tern  becoaes; 

+&**.$&  *  (ci2) 

<5x  dX  w*  wX 

He  neglect  high  order  teras  created  by  the  inner  products  of 
the  eleaents  of  Az<*>,  such  as  the  last  tera  in  relation 
C 12.  The  first  tera  would  be  a  part  of  the  function  P(z<l>) 
in  the  right  side  of  Sguation  CIO  ,  and  the  reaaining  two 
teras  would  be  included  as  part  of  the  expression: 

A  (z<»  >)•  Az<*  >. 


101 


APPENDIX  D 


DESCRIPTION  OP  COMPUTER  PH OGHAMS 


This  appendix  serves  as  a  user  oriented  guide  to  the 
computer  progress  listed  in  App.  E.  The  prograss  are  listed 
in  the  order  used  to  generate  solutions  to  the  problen  for¬ 
mulated  in  this  report.  The  output  consists  of  printed 
results  for  each  iteration  and  punched  output  to  provide  a 
peraanent  file  copy,  which  is  to  be  used  as  input  to  the 
three  plotting  programs  also  listed  in  App.  E.  These  plot¬ 
ting  programs  use  the  Versatec  plotter  and  associated  soft¬ 
ware  *  and  represent  the  only  facility  not  commonly  available 
elsewhere.  The  main  program  SHTfl  as  presented*  is  config¬ 
ured  to  also  produce  output  at  a  monitor  terminal  as  used  on 
CP/CSS  (IBM  360)  or  MVS  (IBM  370).  The  output  in  summary  is 


via  three  devices; 


DEVICE  4  -  typically  a  CRT 

DEVICE  6  -  typically  a  printer 

DEVICE  7  -  typically  a  punch  or  disk  output 


These  can  be  manipulated  in  convenient  configurations  to 


give  a  full  or  partial  (sampled)  display  of  the  solution  at 


each  iteration*  or  merely  to  display  convergence  information 


in  a  simple  format. 
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The  SflTH  prograa  is  the  Bain  prograa  which  together  with 
its  9  subroutines  and  1  IHSL  subroutine  solves  the  non-di- 
aensionalized  fores  of  Egs.  1  through  4.  The  input  to  the 
prograa  consists  of  control  cards  and  data  cards,  naaely; 

1  control  card 

3  or  104  initial  solution  data  cards 
1  control  card 

The  control  card  is  described  here  in  terns  of  the  Fortran 


variable  nanes  thereon. 


IIIII  -  is  the  iteration  counter  (ITBB)  starting  value. 

If  ITBIM  *  1,  then  the  prograa  v3.ll  read  in  the 
next  3  cards  as  data  which  contain  the  short  fora 
starter  solutions.  If  ITBip  is  greater  than  1. 
then  the  prograa  will  read  in  the  next  104  cards  as 
data,  representing  a  detailed  starting  solution. 
This  detailed  starting  solution  is  typically  the 
output  of  a  previous  run  of  the  prograa.  If  ITHIH 
-  00,  then  the  prograa  will  not  read  in  any  data, 
being  thus  inforaed  that  the  problea  is  over. 

ITBAX  -  is  the  last  iteration  value.  The  orograa  ceases 
iterations  at  this  point  and  outputs  the  solutions 
and  auxiliary  inforaation. 

ITHDT  -  enables  the  eaployaent  of  the  net  production  tern 
(BTSXDB)  when  ITBB  is  greater  than  ITHDT,  otherwise 
the  net  production  is  set  to  zero. 

ITBBP  -  for  ITBB  less  than  ITEHP,  the  electron  teapera- 
ture  is  isotheraal,  otherwise  it  is  a  function  of 
E/n.  (see  SHT&  stateaents  899-990)  . 


ITBC 


-  for  ITBB  greater, than  5  and  less  than  ITBC,  the 
wall  boundary  conditions  for  the  species  is 
floated. 


LADJ  -  1  electron  teaperature  is  calculated 
-  0  electron  teaperature  is  isotheraal. 

LP0I  -  1  output  is  punched  on  device  7.  -  9  output  is 

not  punched. 

B  -  sets  the  aesh  spacing  egual  to  H*LA8S  (sheath 
length) . 

FKTB  -is  the  adjustaent  factor  to  the  teaperature  equa¬ 
tion  Xs**  Sec.  17)  and  is  used  when  LADJ  =  1  and 
ITBB  is  less  than  ITBBP.  The  adjustaent  factor  is 
not  used  when  LADJ  is  set  to  zero.  Bhen  the  teap- 
erature  is  coupled  and  LADJ  »  1 ,  then  the  prograa 
calculates  its  own  adjustaent  factor. 

BABB  -  is  reserved  for  user  labelling  and  identification 
which  is  used  at  output. 
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The  solution  set  is; 

PB  voltage  potential  (volts)  . 

SB  electron  density  (10*  •  a-*) 

SI  ion  density  ( 10 4 •  a~3) 

TH  nondisensional  electron  tesperature 
The  short  fora  starter  solution  is  illustrated  in  App. 

E.  It  consists  of  PB,  SB  and  SI  values  at  the  anode,  wall 
and  eguilibriua  boundaries.  The  long  fora  starting  solution 
consists  of  PB,  BE,  SI,  and  TB  distributions,  each  one 
represented  by  25  data  cards  and  a  label  card.  Each  card 
has  7  values,  (see  format-803). 

A  typical  iteration  series  would  start  out  for  about  ten 
iterations  froa  the  short  fora  data  with  ten  volts  at  the 
anode  and  other  appropriate  boundary  conditions.  The  ini¬ 
tial  solution  converges  very  well  for  B  *  0.5,  the  net  pro¬ 
duction  tern  coupled  after  2  or  3  iterations,  and  the  wall 
conditions  not  floated.  This  output  is  then  used  as  input. 
The  voltage  is  anplified  to  say  25  or  35  volts  and  a  new 
control  card  is  structured  to  couple  the  tenperature  or  wha¬ 
tever  else  is  desired.  The  B  value  nay  also  be  changed  to 
reflect  a  saaller  aesh  spacing.  It  aust  be  kept  in  Bind 
that  too  aany  changes  in  one  solution  run  can  be  catas¬ 
trophic. 
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Each  iteration  of  the  solution  by  the  SHTH  program 
accomplishes  the  following; 

(1)  floats  the  HB,  HI  wall  values 

(2)  sweeps  through  the  A-aatrir  and  C-vector  setting 
up  the  equations  that  satisfy  the  boundary  condi¬ 
tions  and  the  Internal  mesh  conditions 

(3)  couples  the  electron  temperature  equation 

(4)  couples  the  net  production  term 

(5)  calculates  the  anode  electron  density  required  to 
satisfy  current  matching  btween  the  anode  and  the 
f ree-streaa. 

The  system  of  equations  has  the  fora  presented  in  Eq. 

CIO  in  App.  C,  and  is  solved  by  use  of  the  IHSL  routine 
LEQT1B.  The  solution  vector  is  then  updated  from  the  return 
information  contained  in  vector  c. 

There  are  two  solution  convergence  monitors.  The  first 
is  the  sum  of  the  squares  of  the  differences  (spectral  norm) 
of  the  -P(z<l>)  values  from  Eg.  CIO.  The  other  is  the  spec¬ 
tral  norm  of  the  Az(t>  upon  return  from  LEQT1B.  Both  of 
these  are  available  at  the  video  terminal  (device  4)  and 
printer  (device  6). 

After  solution  update,  the  elctron  temperature  is  calcu¬ 
lated  and  adjusted  per  the  control  card  options.  It  should 


be  pointed  oat  that  on  the  last  iteration,  the  teaperature 
distribution  is  autoaaticaily  calculated  for  output,  but  is 
not  used  in  the  scheme,  unless  otherwise  coupled  per  the 
control  card. 

The  matching  of  the  anode  current  to  the  free-stream 
current  is  accomplished  by  solving  for  the  anode  electron 
density  using  subroutine  MODE.  The  call  to  this  subprogram 
will  require  a  change  to  the  SHTH  program  as  it  is  presently 
overridden  with  a  comment  card,  (prior  to  statement  9000) . 

The  final  outputs  are  straight-forward,  subroutines 
JWALL  and  JOULE  are  called  in  this  phase  of  the  program,  in 
order  to  compute  the  current  into  the  wall  and  the  Joule 
heat  distribution.  The  electric  field,  net  production  and 
space  charge  density  distributions  are  also  calculated  and 
provided  as  punched  and  printed  output. 

The  punched  output  is  used  as  input  for  the  STBEAH,  PICT 
and  PLOT  programs.  These  programs  use  the  Tersatec  plotter 
to  provide  a  presentation  of  the  current  stream  lines,  an 
oblique  perspectives  of  the  solutions  and  plots  of  the  solu¬ 
tions  as  they  vary  normal  to  the  wall  extending  from  the 
anode  site.  The  details  of  these  programs  are  combersome, 
and  it  is  wise  to  leave  them  to  the  user  for  inspection. 
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The  first  two  distributions  in  the  punched  output  are 
used  in  the  STBB&M  prograa,  requiring  that  the  nonpertinent 
data  be  resowed.  The  next  8  distributions  are  used  for  the 
PLOT  and  PICT  programs.  The  PICT  prograa  has  an  option  card 
filled  with  ones  or  zeros  in  the  first  eight  positions. 

This  option  card  enables  the  user  to  elect  which  distribu¬ 
tions  to  plot.  The  plots  provided  in  this  report  as  Figs. 

12  through  14,  were  produced  using  the  PICT  and  PLOT  pro¬ 
gress.  The  progress  are  sisple  and  speak  for  thesselves. 


APPENDIX  E 
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COMPOTES  PROGRAM  LISTINGS 

This  appendix  contains  a  listing  of  the  programs  used  in 
the  solution  of  the  equations  presented  in  this  work 
(namely:  Bqs.  9,10  and  11),  and  the  auxiliary  programs  that 
provided  the  plots  for  the  report. 


SHTH  (MAIN  PROGRAM) -  109 

DELOP -  127 

GRDNT -  129 

TEMP -  131 

RTSIDE -  132 

PATCH -  133 

ANODE -  134 

JBALL -  136 

JOOLE -  137 

OOTPOT - — -  138 

PLOT -  139 

PICT -  141 

STREAM -  144 

INITIAL  INPOT  (sample) -  150 


s 
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